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Abstract - We report a multiple-site mean-field 
analysis of the zero-temperature phase diagram 
for ultracold bosons in realistic optical superlat- 
tices. The system of interacting bosons is de- 
scribed by a Bose-Hubbard model whose site- 
dependent parameters reflect the nontrivial pe- 
riodicity of the optical superlattice. An analytic 
approach is formulated based on the analysis of 
the stability of a fixed-point of the map defined 
by the self-consistency condition inherent in the 
mean-field approximation. The experimentally 
relevant case of the period-2 one-dimensional su- 
perlattice is briefly discussed. In particular, it is 
shown that, for a special choice of the super- 
lattice parameters, the half-filling insulator do- 
main features an unusual loophole shape that the 
single-site mean-field approach fails to capture. 
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1. INTRODUCTION 

Although originally introduced for liquid He- 
lium in confined geometries pQ, the Bose- 
Hubbard (BH) model proves successful in de- 
scribing ultracold atoms trapped in optical lat- 
tices [2]. In this framework, the sites of the am- 
bient lattice correspond to the local minima of 
the effective optical potential created by counter- 
propagating laser beams, and the height of the 
potential barriers between adjacent minima, pro- 
portional to the laser intensity, determines the 
hopping amplitude in the BH model. Such a 
direct relation allows an unprecedented exper- 
imental control of the model parameters, and 
plays a key role in experiments aimed at reveal- 
ing the superfluid-insulator transition character- 
izing the BH model 

In general, the superposition of simple opti- 
cal lattices with commensurate lattice constants 
gives rise to confining potentials characterized 
by a richer periodicity, the so-called superlattices 
Reference JSJ reports on a recent experiment 
where a simple ID optical superlattice is cre- 
ated by superimposing two ID optical lattices 
j^j and a cigar-shaped magnetic potential pro- 
viding confinement in the transverse direction. 
Each of these optical lattices is created by the 
interference pattern of two laser beams crossing 
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at a given angle. The lattice constants, deter- 
mined by such crossing angle, is chosen to be d\ 
and d,2 = 3di, so that the supercell of the result- 
ing optical potential contains three local minima. 
Following this scheme, a ^-periodic ID superlat- 
tice — i.e. a lattice characterized by a £-site 
supercell — can be created by a suitable adjust- 
ment of the crossing angles determining d\ and 
d 2 . 

Following the tight-binding-like approach of 
Ref. |2j a system of ultracold alkali atoms con- 
fined in a ID optical superlattice comprising M 
sites is described by the following BH Hamilto- 



man, 



M jj 

H = [ -77 n k(n k - 1) - (fj, - v k )n k 



k=l 



tk{aka k+1 + a k dk+i] 



(1) 



where a£, a k and n k = a^a k are respectively 
the boson creation, annihilation and number op- 
erators relevant to the site labeled k. As to the 
Hamiltonian parameters, U > accounts for on- 
site repulsion (proportional to the atomic scat- 
tering length) , // is the grand canonical chemical 
potential, v k is the local potential at site k and t k 
is the hopping amplitude between adjacent sites 
k and k + 1. The i-site periodicity of the super- 
lattice yields 



tk+sl = tk = tT k , Vk+ S t = Vk = Wk, 



(2) 



where s = 0, . . . ,M/£ — 1 labels the supercells 
and t, v are scaling coefficients directly related 
to the intensity of the laser beams giving rise to 
the optical potential. 

As it is well known, in the homogeneous case 
£ = 1, Hamiltonian |JTJ) is characterized by 
the superfluid-insulator quantum phase transi- 
tion pQ. More in detail, the competition be- 
tween the on-site repulsion and the kinetic en- 
ergy — proportional to U and t, respectively — 
gives rise to a zero-temperature phase diagram 
in the (i/U-t/U plane consisting of an extended 



superfluid phase and a series of adjacent Mott- 
insulator lobes. In the latter the system is re- 
markably characterized by a commensurate pop- 
ulation, i.e. by an integer filling. Several numeri- 
cal and analytical approaches have been adopted 
for the study of such zero-temperature phase di- 
agram. We refer the reader to Ref. [7j for a brief 
review of such techniques. 

Recently, some attention has been devoted to 
the phase-diagram of superlattice BH models 
(HI H3 El ITT] - In general, incompressible Mott 
domains are expected to occur in correspondence 
of critical fractional filling. In the case of ID £- 
periodic superlattices such critical fillings are in- 
teger multiples of £ _1 . Furthermore it has been 
shown that, when the local potentials Vj are not 
all different from each other, some of the Mott 
domains exhibit an unusual loophole shape 

In this paper we study the zero-temperature 
phase-diagram of Hamiltonian (JTJ adopting a 
multiple-site mean-field approach generalizing 
the technique introduced in Ref. ■ The latter 
provides satisfactory qualitative results for the 
quantum phase transition occurring in the ho- 
mogeneous case, but fails to predict the loophole 
insulator domains that may appear in the case 
of superlattices . We mention that a two-site 
mean-field approach is adopted in Ref. [7] for 
the study of homogeneous lattices. We further- 
more show that, in general, the zero-temperature 
phase diagram can be worked out analyzing the 
stability of a particular fixed-point of the map 
defined by the mean-field self-consistency condi- 
tion ^Oj. Most of such analysis can be carried 
out analytically based on a perturbative expan- 
sion of the spectrum of the mean-field Hamil- 
tonian. This allows to determine the phase di- 
agram by solving a numerical problem that is 
much less demanding than the iterative proce- 
dure standardly used to deal with the original 
self-consistency equations. Furthermore, in some 
special cases, entirely analytical results can be 
obtained. In particular, our method provides the 
analytical description of the Mott-lobe bound- 
aries of the homogeneous case E] • Exploit- 
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ing our method we analyze the realistic case of 
a t = 2 ID superlattice. In particular, we study 
the insulator domain corresponding to the crit- 
ical filling / = 1/2, showing that its usual lobe 
shape shrinks at the bottom and turns into a 
loophole as the potential offset vi — v\ between 
the sites of the supercell vanishes. Furthermore, 
we provide the exact analytic description of the 
boundaries of such loophole domain. 

2. MULTIPLE-SITE MEAN-FIELD 

In the simple case of the usual lattice, 
I = 1, qualitative information about the zero- 
temperature phase diagram of Hamiltonian Q 
can be obtained making use of the single-site 
mean-field approach introduced in Ref. |13j . De- 
noting (•) the expectation value on the ground 
state, it is assumed that for every k, 

a kat +l = {a k )al +l + a k {al +l ) - (a k )(a+ +1 ). (3) 

This allows to recast Hamiltonian (0) as the sum 
of M single-site Hamiltonians, H J2k=i^k, 
where 



U 



nk(n k - 1) - iin k 



t(a k+1 + a k -i)(a k + - a k ) (4) 



and the so-called superfluid parameters are to be 
determined self-consistently as 1 



a k = (a k ) = (a£). 



(5) 



After the translational invariance characterizing 
the system is taken into account, a k = a k +i = 
a, the Hamiltonians in @ are decoupled and 
become formally identical. The original prob- 
lem hence reduces to the study of one single- 
site Hamiltonian. In this framework, the Mott- 
insulator domains in the /i/U — t/U phase dia- 
gram are characterized by the vanishing of the 
superfluid order parameter. Indeed in this case it 



1 The superfluid parameters defined in Eq. JSJ are real 
since the boson operators in Eq. |g) have a real represen- 
tation on the usual Fock basis. 



is easy to check that local density of bosons (n k ) 
is pinned at an integer value and the system is 
incompressible. 

On inhomogeneous structures translational in- 
variance is lost, and the single-site Hamiltonians 
in Eq. ((3} are coupled by the self-consistency 
conditions This single-site mean-field ap- 

proach gives fairly satisfactory qualitative results 
for superlattices whose supercell features local 
potentials v k all different from each other |1U| . 
but fails to capture the loophole-shaped insula- 
tor domains appearing when some of these po- 
tentials are equal |12| . 

A more structured approach, taking into ac- 
count the non-trival periodicity of a ^-periodic 
superlattice, consists in adopting approximation 
© every ^-th site. By so doing, Hamiltonian P| 
becomes the sum of identical i-site Hamiltoni- 
ans, one for each supercell, and, as in the single- 
site approach, the original problem reduces to 
the study of one of such supercell Hamiltonians. 
Dropping the supercell index, the latter reads 



H 



T^n k (n k -!)-(//- v k )n k 



E 

fc=i 

e-i 

E t k( a ka k+ i + a k a k+ i) 
k=l 

t£ |a^(ai + af) + a\(a£ + a^) — 2a\ai 

(6) 



where 



aj = (dj) = {dj} + , j = 1, 



(7) 



This approach is expected to give satisfactory re- 
sults if approximation (j3J) is adopted for the hop- 
ping terms characterized by the lowest hopping 
amplitude, ti < t^- As in the single-site approx- 
imation, the Mott-insulator phase is character- 
ized by vanishing superfluid parameters, ot\ = 
a£ = 0. In this situation, the mean-field Hamil- 
tonian 7i commutes with the total number of 
bosons (in the supercell), J2 k =i n k- Hence the 
expectation value of the latter on the ground- 
state is fixed to an integer value determined by 
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the Hamiltonian parameters and, quite interest- 
ingly, the filling of the system / = M~ l J2h=i n k 
is a multiple of l^ 1 . 

The most standard approach to the mean-field 
problem defined by Eqs. © and (JJJ) consists of 
an iterative numerical procedure. More in de- 
tail, the superfluid parameters appearing in the 
multiple-site Hamiltonian @ at a given iteration 
are determined evaluating Eqs. (J7J) on the ground 
state of the previous iteration. The procedure is 
arrested when the value of the superfluid param- 
eters does not change significantly between two 
subsequent iterations. 

3. ANALYTICAL APPROACH 

The standard iterative procedure illustrated 
above shows that solving the self-consistency 
problem in Eqs. (JSJ) and (J7J) amounts to finding 
a stable fixed point of the map 



a'x = Fi(ai,ae) 



(8) 



where Fj(a\,ai) = (aj), j = 1,1. Note that the 
choice a.\ = ct£ = 0, corresponding to the Mott- 
insulator phase, is a fixed point of the map in 
Eq. (JHJ) for any value of the Hamiltonian param- 
eters fi, U, {tfc}, {ffc}- Indeed, as we mention in 
the previous section, in this situation the ground 
state of the system belongs to a fixed-number 
subspace, and the expectation values in Eq. (JJJ) 
necessarily vanish. This means that the insula- 
tor domains are characterized by choices of the 
Hamiltonian parameters making the fixed point 
a% = a.g = stable. According to the standard 
criterion, this happens when the absolute value 
of the eigenvalues of the Hessian matrix 



d ai Fi(ai,a e ) 



d at Fi(ai,ae) 



d ai F e (ai,a e ) d ae F e (a>i,> 



0) 



is smaller than 1. 

Note that the Hessian matrix is completely de- 
termined by a first order expansion of Fj in the 
parameters ct\, ag, which can be in turn obtained 



from a first order expansion to the ground state 
of TC in the same parameters. Since the term 
tta,\ai in Eq. © does not contribute first order 
corrections to the ground state of Tt, it can be 
discarded without loss of generality. After doing 
this, the desired first order approximation can be 
obtained using tt as the perturbative parameter, 
since it multiplies all of the first order terms in 
a\ and ag appearing in Eq. ©: 



H 
V 



n + tgV 

— Q^(ai + af 



(10) 
(11) 



Since the unperturbed Hamiltonian T~Lq com- 
mutes with the total number of bosons (in the 
supercell), its eigenstates belong to fixed num- 
ber subspaces. Denoting \(f>h) such eigenstates 
and e/j the relevant eigenvalues, the first or- 
der approximation to the ground state of 7i is 
IV') ~ I0o)+^|^/ ; i)) where |</>o) is the unperturbed 
ground-state and 



eo - £h 



(12) 



This means that 



Fj{a>i,a e ) 



{if)\aj\ip) « (^ |%#i) + (ipi\aj\(f)o) 
aiCje + agCji, (13) 



where 



Cjk '- 



') h \aj+a'j\(j) )((j) h \ak+a^ \<p ) 
e/i - e 



(14) 



Therefore, as we mention, the Hessian matrix in 
Eq. Q is determined in terms of the coefficients 
appearing in the first order approximation (|13|) . 
defined in Eq. ()14|) . and the condition for the 
stability of fixed point a\ = at = is 



\cu ± y/ciicee\ < 1 



(15) 



Since the coefficients in Eq. (|14j) depend on the 
Hamiltonian parameters in Eqs. (^Q) and (J5J) 
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through the eigenvalues and eigenstates of Ho, 
inequality (|15j) allows to determine the regions 
of the n/U-t/XJ plane pertaining to the Mott- 
insulator phase. According to the above discus- 
sion, within such phase the (integer) number of 
bosons in each supercell is N = (0o| J2k=i n k\<Po) 
and corresponds to the fractional filling / = N/l. 

Note that the study of the phase space by in- 
equality (fT3|) is much less demanding than the 
standard iterative procedure briefly illustrated in 
the previous section. Indeed, for a given choice of 
the Hamiltonian parameters, the latter involves 
the iterative diagonalization of a matrix whose 
size is Y,k=odk, where d k = is the di- 

mension of the subspace relevant to k bosons in 
£ sites and C provides a cutoff for the in principle 
infinite Hilbert space of the problem. 

Conversely, no iterative procedure is required 
for the study of inequality (|T5|) . Indeed it is suffi- 
cient to diagonalize only the three (independent) 
blocks of Ho relevant to the total numbers of 
bosons N — 1, N and N + 1, where N is the cell 
population characterizing the Mott domain un- 
der investigation. Furthermore, as we illustrate 
in the following, in some simple cases inequal- 
ity (0 can be studied in a completely analytical 
way. 

4. RESULTS: 1 = 2 SUPERLATTICE 

In this section we consider the realistic case of 
a t = 2 ID superlattice, that can be created as 
in Ref. [Hj by superimposing two homogeneous 
ID optical lattices with lattice constants d\ and 
d 2 = 2d\. The insulator domains relevant to the 
lowest fractional fillings (dark gray), as evalu- 
ated by means of a numerical study of inequal- 
ity (fH)|) , are displayed in Figure ^ for the pa- 
rameter choice U = 1.0, t\ = 1.0, t 2 = 0.3, 
vi = 0.00, and, from top to bottom panel, 
v 2 = 0.12, 0.06, 0.03, 0.00. Note that the width 
at t/U = of the half- filling insulator domain 
equals the energy offset between the sites of the 
same supercell, v 2 — v\ ^U]. As the latter van- 
ishes, such insulator domain assumes an unusual 
loophole shape [T2] . 




Figure 1: Superfluid "clamping" around a frac- 
tional insulator domain for a I = 2 superlattice 
BH Hamiltonian. Each panel shows the region 
of the phase diagram containing the insulator 
domains (dark gray areas) relevant to the lower 
critical fillings, also shown. The insets contain 
a pictorial representation of the relevant effec- 
tive optical potential. As the energy offset be- 
tween the lattice sites in the same supercell de- 
creases — from top to bottom — , the superfluid 
phase (white) clamps around the half-filling in- 
sulator domain, which assumes an unusual loop- 
hole shape. 



In this special case, (v\ = v 2 , bottom panel 
of Fig. , the study of inequality (|T5|) involves 
the diagonalization of 2 x 2 matrices, and can 
be carried out in a completely analytical way. 
After some calculations, it is possible to show 
that the loophole domain border is determined 
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by the following equation 

= ?>t\ (n+T 2 )t 3 

+ [C/(2r 2 -n)+M (5r 1 + 2r 2 )]r 1 i 2 
+ /u 2 (n -r 2 ) t + Ufi 2 -n 3 (16) 

A simple analysis shows that the loophole do- 
main disappears when no positive t satisfies the 
preceding equation for \i = 0, i.e. when t\ < 2 r 2 . 
However, this is an artifact introduced by the 
mean-field approximation, which is known to 
provide at best qualitative information. As it is 
shown in Ref. |T2|, the loophole domain can be 
proven to exist for any r 2 ^ T\ by resorting to 
the exact mapping between the hard-core limit 
(t/U <C 1) of Hamiltonian Q and the model for 
spinless non-interacting fermions on the same su- 
per lattice. 

5. CONCLUSIONS 

In this paper we introduce a multiple-site 
mean-field approach to the study of the zero- 
temperature phase diagram for ultracold bosons 
in realistic one-dimensional superlattices. A per- 
turbative expansion in the hopping amplitudes 
between neighbouring supercells allows to recast 
the self-consistency costraints involved in this 
approach into a problem that requires a numer- 
ical effort much less demanding than the usual 
iterative procedure. 

Relying on such a multiple-site mean-field ap- 
proach, we supply some explicit results for the 
experimentally relevant case of a 2-periodic su- 
perlattice. In particular we show that, as the 
energy offset between the sites of the supercell 
decreases, the superfluid phase clamps around 
the half-filling insulator domain, which assumes 
an unusual loophole shape for vanishing energy 
offset. Our multiple-site mean-field approach 
shows that this "clamping effect" of the super- 
fluid phase likewise occurs around all of the frac- 
tional filling insulator domains of the I = 2 
superlattice. Similar, and even more complex 
"clamping effects" can be shown to occur on 
more structured superlattices II T. . 
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